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Abstract 

Peskin’s  Immersed  Boundary  (IB)  method  is  one  of  the  most  popular  numerical  methods 
for  many  years  and  has  been  applied  to  problems  in  mathematical  biology,  fluid  mechanics, 
material  sciences,  and  many  other  areas.  Peskin’s  IB  method  is  associated  with  discrete  delta 
functions.  It  is  believed  that  the  IB  method  is  first  order  accurate  in  the  L°°  norm.  But  almost 
none  rigorous  proof  could  be  found  in  the  literature  until  recently  [9]  in  which  the  author  showed 
that  the  velocity  is  indeed  first  order  accurate  for  the  Stokes  equations  with  a  periodic  boundary 
condition.  In  this  paper,  we  show  the  first  order  convergence  with  a  log  h  factor  of  the  IB  method 
for  elliptic  interface  problems  essential  without  the  boundary  condition  restrictions.  The  results 
should  be  applicable  to  the  IB  method  for  many  different  situations  involving  elliptic  solvers  for 
Stokes  and  Navier-Stokes  equations. 
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discrete  Green  function,  discrete  Green’s  formula. 
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1  Introduction 

Since  its  invention  in  1970’s,  the  Immersed  Boundary  (IB)  method  [10]  has  been  applied  to  almost 
everywhere  in  mathematics,  engineering,  biology,  fluid  mechanics,  and  many  many  more  areas,  see 
for  example,  [11]  for  a  review  and  references  therein.  The  IB  method  is  not  only  a  mathematical 
modeling  tool  in  which  a  complicated  boundary  condition  can  be  treated  as  a  source  distribution 
but  also  a  numerical  method  in  which  a  discrete  delta  function  is  used.  The  IB  method  is  robust, 
simple,  and  has  been  applied  to  many  problems. 
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It  is  widely  believed  that  Peskin’s  IB  method  is  only  first  order  accurate  in  the  L°°  norm. 
However,  there  was  almost  none  rigorous  proof  in  the  literature  until  recently  [9],  in  which  the 
author  has  proved  the  first  order  accuracy  of  the  IB  method  for  the  Stokes  equations  with  a  periodic 
boundary  condition,  see  also  [3]  in  which  the  author  showed  that  the  pressure  obtained  from  IB 
method  has  0( h 1^2)  order  of  convergence  in  the  L 2  norm  for  a  ID  model.  However,  there  are  few 
theoretical  proof  on  the  IB  method  for  elliptic  interface  problems.  This  is  the  main  motivation  of 
this  section.  The  main  goal  of  this  paper  is  to  provide  a  convergence  proof  for  the  IB  method.  We 
will  show  that  with  reasonable  assumptions  (zero-th  moment  condition  and  first  order  interpolation 
property)  on  the  discrete  delta  function  used  in  the  IB  method,  the  IB  method  is  indeed  first  order 
convergent  in  the  L°°  norm.  Our  proof  is  essentially  independent  of  the  boundary  conditions  and 
it  is  valid  in  ID,  2D,  and  3D  cases.  The  result  should  be  applicable  for  many  IB  methods  involving 
Stokes  and  Navier-Stokes  solvers. 


2  Proof  of  the  convergence  of  the  IB  method  in  ID 


We  will  give  a  proof  for  the  ID  model  first, 

u"  =  cS(x  —  a)  0  <  x  <  1,  0  < 

in  this  section.  Note  that  the  analytic  solution  to  the 

—cx  (1  —  a) 


i(x)  = 


— ca  (1  —  x) 


a  <  1,  it(0)  =  u(  1)  =  0, 

equation  above  is 

if  x  <  Xk, 

otherwise. 


(1) 

(2) 


Given  a  uniform  Cartesian  grid  x*  =  ih,  i  =  0, 1,  •  •  •  ,  n,  h  =  1/n,  the  IB  method  leads  to  the 
following  system  of  linear  equations, 

Ui-i  —  2 Ui  +  Ui- |-i  x  !  \  ■  i  o  i  ro\ 

- — - =  cdh{Xi-a),  i  =  1,2,-  •  •  ,n  -  1,  (3) 

where  Ui  is  the  finite  difference  approximation  of  the  solution  it(xj),  and  5h{xi  —  a)  is  a  discrete  delta 
function  applied  to  the  grid  point  x*.  In  the  matrix- vector  form,  the  above  finite  difference  equations 
can  be  written  as  ATJ  =  F,  where  A  is  the  tri-diagonal  matrix  formed  by  the  discrete  Laplacian.  It 
is  well  known  that  —A  is  a  symmetric  positive  definite  matrix  and  diagonally  dominant.  Note  that, 
a  discrete  delta  function  has  a  compact  support  in  the  neighborhood  of  interface,  that  is,  Sh(x)  ^  0 
only  if  |x|  <  Wh,  where  W  is  a  constant.  Commonly  used  discrete  delta  functions  include  the  hat 
discrete  delta  function  ( Shat(x )  with  W  =  1)  and  the  cosine  discrete  delta  functions  (Scosme(x)  with 
W  =  2),  see  for  example,  [8].  Note  that,  when  we  use  the  hat  delta  function,  the  result  is  the  same 
as  that  of  the  IIM  for  the  simple  model.  The  solution  to  the  finite  difference  equations  is  the  same 
as  true  solution  if  there  are  no  round-off  errors,  that  is 


u{xi- 1)  -  2 u(xi)  +  u(xi+ 1) 
h2 


=  c6hhat  (Xi 


a) , 


*  =  1,2,---  ,  n  —  1, 


(4) 
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see  for  example,  [1,  8].  But  this  is  not  the  case  for  other  discrete  delta  functions. 

We  define  the  error  vector  as  E  =  {E,},  where  Ei  =  u(xi )  —  t/j.  The  local  truncation  error  is 
defined  as  T  =  {T(}, 


Ti  = 


i{xi- 1)  -  2 u{xi)  +  u(xi+ 1) 
h? 


-  c6h  (xi  -  a) . 


(5) 


With  the  definition,  we  have  =  F,  Ah u  =  F  +  T,  and  therefore  Ah  E  =  T.  For  the  hat  discrete 
delta  function,  we  have  |T*|  =  0  for  all  Vs.  For  other  discrete  delta  functions,  we  have  \Tj\  =  0(l/h) 
for  a  few  grid  points  neighboring  the  interface  a.  So  the  interesting  question  is:  Why  is  the  IB 
method  still  first  order  accurate,  that  is,  H-Elloo  =  0(/i)?  To  answer  this  question,  we  first  introduce 
the  following  lemma. 

Lemma  2.1.  Let  Ahy  =  ek,  where  ek  is  the  k-th  unit  base  vector,  yo  =  yn  =  0,  then 

(  - hxi  (1  -  xk)  if  i  <  k, 

Ui  =  \  (6) 

[  — hxk  (1  —  Xi)  otherwise. 


The  significance  of  this  lemma  is  that  the  solution  is  order  h  smaller  than  the  concentrated 
source. 

Proof:  We  note  the  following  identity 

Ahy  =  h^  =  hdhhatek.  (7) 

n 

For  this  simple  case,  the  IB  and  IIM  are  identical.  Thus  from  the  Immersed  Interface  Method,  see 
[1,  8],  we  know  that  y  is  the  exact  discrete  solution  at  the  grid  points  of  the  following  boundary 
value  problem 

•u(O)  =  u(  1)  =  0,  (8) 

if  x  <  xk , 

t  .  (9) 

otherwise. 

This  completes  the  proof.  Note  that  \y(x)\  <  h.  From  this  lemma,  we  have  the  following  corollary. 

Corollary  2.2.  Let  Ahy  =  r,  then  |j/(a:)|  <  hW,  where  W  is  the  number  of  non- zero  components 
of  r. 

Notice  that  for  a  discrete  delta  function,  it  should  satisfy  at  least  the  zeroth  moment  equation, 
see  [2],  that  is 

y ^Sh(xj  -  a)  =  1,  (10) 

i 

corresponding  to  the  continuous  case  f  Sh{x  —  a)dx  =  1.  Now  we  are  ready  to  prove  the  main 
theorem. 


whose  solution  is 


u"  =  h5(x  —  Xk),  0  <  x  <  1, 
[  -hx(l-xk) 

y(x )  =  < 

I  -hxk  (1  -  x) 
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Theorem  2.3.  The  Immersed  Boundary  method  for  (1)  is  first  order  accurate,  that  is 


llEHoo  <Ch, 


(11) 


where  C  is  a  constant. 


Proof:  We  can  decompose  the  local  truncation  error  into  two  groups 

T  =  Tre®  +  T  irreB  (12) 

where  llT^Hoc  =  0  corresponds  to  the  local  truncation  errors  at  regular  grid  points  where  Sh(xi  — 
a)  =  0  and  the  true  solution  is  piecewise  linear.  Note  that,  we  have 

n— 1 

Ti  =  Treg  +  E  rj\irreg  =  Q  +  E  Tirreg.  (13) 

On  the  other  hand,  we  also  have 

n—  1  n—  1  n—  1 

c5%at(xi  -  a)  -  cSh(xi  —  oi)  =  0,  (14) 

2=  1  2=1  2=1 

since  the  finite  difference  method  using  the  discrete  delta  function  gives  the  exact  solution  at  the 
grid  points.  Thus  we  have  =  T[rreg  =  0.  We  divide  T-rreg  into  two  groups,  one  with  all 
positive  Tf s  denoted  as  j'*rre9<+-  the  other  one  is  all  the  negatives  denoted  as  X-,reg’^.  Since 
Tf>reg,+  _|_  y^T *rre9-  _  rpirreg,+  an(j  rpvrreg-  mugt  have  the  same  order  of  the  magnitude  ( 

0(l/h))  although  those  index  i  are  different  except  that  | —  a\  <  Wh  is  true  for  all  irregular  grid 
points.  Because  the  solution  is  linear  with  c,  we  have 

E  =  Af1  T  =  Af1  (Tirres’+  +  Tirreg’~)  .  (15) 

From  the  solution  expression,  we  know  that,  assuming  that  Xk  <  a, 


Ek 


- xk 


T-  (1 

i  3 


=  -xk{  1  -  a) 
=  0(Wh), 


^  rjiirreg,+  _|_  ^  ^  rpirreg , 


+  0(Wh) 


after  we  expand  all  xfs  and  Xk  s  at  a  and  since  all  related  Xi  and  Xj  are  within  Wh  distance  from 
the  interface  a.  The  proof  when  Xk  >  a  is  similar  except  that  we  need  to  use  the  solution  for  x  >  a. 
This  completes  the  proof. 
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3  Proof  of  the  convergence  of  the  IB  method  in  2D 


The  discussion  for  2D  problems  is  much  more  challenging  since  the  interface  is  often  a  curve  instead 
of  a  point.  In  [9],  the  author  has  proved  the  first  order  accuracy  of  the  IB  method  for  the  Stokes 
equations  with  a  periodic  boundary  condition  in  2D.  However,  there  are  almost  no  theoretical  proof 
on  the  IB  method  for  elliptic  interface  problems  or  other  general  boundary  conditions.  We  will  prove 
that  the  result  obtained  from  the  IB  for  the  elliptic  interface  problem  is  indeed  first  order  accurate 
in  2D  as  well  in  this  section. 


Consider  the  2D  elliptic  interface  problem 


A u(x,y)  =  f{x,y)  +  fr  v(s)6  (x-X(s))  ( y-Y(s))ds 
u(x,y)\dQ  =u0(x,y), 


(16) 


where  we  assume  that  /  £  C(fl),  T  £  C,  v(s)  £  C.  Without  loss  of  generality,  we  assume  that  H 
is  a  unite  square  0  <  x,y  <  1.  The  problem  can  be  decomposed  as  the  sum  of  the  solutions  of  the 
following  two  problems.  The  first  one  is 


Aui  (x,y)  =  f(x,y ) 
ui(x>y)\dti  =  uo  (x>y)> 

which  is  a  regular  problem  whose  solution  u\(x,  y)]  £  C2(f2).  The  second  problem  is 


(17) 


A u2(x,y)  =  J  v(s)S(x-X(s))(y-Y(s))ds 
Mx,y)\an  =0, 

The  solution  to  the  second  problem  is  equivalent  to  the  following  problem 

~du2 1 


Au2(x,y)=0,  [u2]r  =  0, 


dn 


=  v(s), 


(18) 


(19) 


U2(x,y)\dn  =  0. 

The  solution  to  the  original  problem  is  u  =  u\  +  u2.  Since  u i  is  the  solution  to  a  regular  problem, 
it  is  enough  just  to  consider  u2(x,y),  and  we  will  simply  use  the  notation  u(x,y)  for  u2(x,y). 

The  Peskin’s  IB  method  for  the  problem  including  the  following  steps. 


•  Generating  a  uniform  Cartesian  mesh  Xi  =  ih ,  yj  =  jh ,  i,  j  =  0, 1,  •  •  •  ,  n. 

•  Replace  the  partial  derivatives  with  the  finite  difference  approximation  and  use  a  discrete  delta 
function 


Ui-i,j  +  Ui+ ij  +  Uitj- 1  +  Uij. (_i  —  4 Uij  _  r<IB 

h 2  _  ij 


i,j  =  0, 1,  -  -  -  ,n, 


(20) 


Nb 

C-f  =  ^2  v^h  (Xi  -  xk)  Sh  {yj  -  Yk)  A Sfc, 

fc= l 


(21) 
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where  (Xk,  Yk),  k  =  1,  2,  •  •  •  ,  7V&  is  a  discretization  of  the  interface  T,  and  Vk  ~  v(sk). 

•  Solve  the  finite  difference  system  of  equations  above  to  get  an  approximation  solution  {  Ul3  } . 


3.1  Assumptions  and  the  discrete  delta  function 


To  guarantee  enough  accuracy,  we  assume  that  maxfcjsfc}  =  As  ~  0(h ),  and 

n—  1  Nb  p 

y.  h 2  y  vkSh  (Xi  -  Xk)  Sh  (: Uj  -  Yk)  A sk  =  /  v(s)ds  +  0(h), 

i,j=l  fc= 1 

see  for  example,  [9],  which  corresponds  to 


n  Jr 


)(s)5  (x  —  X(s))  5  (y  —  Y(s))  dsdxdy  =  J  v(s)ds. 


(22) 


(23) 


From  f fn  u(x,  y)S  (x  —  X)  6  (y  —  Y)  dxdy  =  u(X,  Y),  we  should  also  have  the  interpolation  property, 


n—  1 


y  h2u(xi,yj)Sh  (xi  -  Xk)Sh  (■ yj  -  Yk)  =  u(X,Y)  +  O(h). 
i,j= i 


(24) 


The  discrete  delta  function  has  a  compact  support,  that  is, 

5h(xi-Xk)  =  0,  if  \Xi-Xk\>Wh,  and  Sh(yj-Yk)  =  0,  if  \Vj  -  Yk\  >  Wh,  (25) 
where  x.y  =  (xi,yj)  and  IT  is  a  constant. 

We  define  the  error  vector  as  E  =  {Eij},  where  Ei3  =  u(xi,yj)  —  Uij.  The  local  truncation  error 
is  defined  as  T  =  {TV,}, 

u(xi-!,yj)  +  u(xi+ 1),  yj)  +  u{xu  Uj- 1)  +  ufa,  yj+ 1)  -  4  u[xi,  yj)  IB 


T  — 


h 2 


° ij  ‘ 


(26) 


In  the  matrix  vector  form,  we  have  =  F,  A/>u  =  F  +  T,  and  therefore  y4/,E  =  T,  where  Ak 

is  the  matrix  formed  by  the  discrete  Laplacian.  In  general,  we  have  | Ttl  \  =  0(1  /h)  for  grid  points 
neighboring  the  interface  T  except  for  the  correction  terms  using  the  Immersed  Interface  Method 
(IIM)  [6,  7,  8]  for  which  we  have  \TEm\  =  0(h)  at  irregular  grid  points  where  the  interface  cuts 
through  the  standard  five-point  stencil. 

So  it  is  interesting  that  the  local  truncation  errors  can  have  order  0(1 /h)  at  some  grid  points, 
but  the  global  error  is  still  of  order  0(h).  There  has  to  be  some  kind  of  cancelations  of  the  errors, 
which  can  be  seen  from  our  proof  process. 

Definition  3.1.  Let  etJ  be  the  unit  grid  function  whose  values  are  zero  at  all  grid  points  except  at 
x.y  =  (a ’i,yj),  =  1.  The  discrete  Green  function  value  at  x/m  centered  at  Xy  is  defined  as  the 
solution  of  the  following 


G  (Xij,X/m)  —  (  (Ah)  Gij  2 


(27) 
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See  for  example,  [1,  4,  5,  9]  for  more  discussions  about  the  discrete  Green’s  functions.  To  prepare 
the  convergence  proof,  we  first  list  some  lemmas  that  are  used  in  the  proof. 

Lemma  3.2.  The  discrete  Green’s  formula  and  an  error  estimate.  Assuming  that  the  distance 
between  T  and  <9f l  is  0(1),  that  is,  dist(T,i 9fi)  ~  0(1),  then  we  have 

'V'  A hu(xi,yj)h2  =  [  Tr-ds  +  0(h)  =  [  v{s)ds  +  0(h),  (28) 

Jan  on  Jr 


where 


A  hu{xi,yj)  = 


+  u(xi+i),yj)  +  u(xi,yj-i)  +  u(xi,yj+1)  -  4 u(xi,yj) 
h 2 


(29) 


is  the  discrete  Laplacian  using  the  standard  five-point  stencil. 


Proof:  We  first  show  the  discrete  Green’s  formula  by  expanding  the  summation.  After  cancela¬ 
tion  of  interior  terms,  only  boundary  terms  are  left.  Thus,  we  get 


E  A  hu(xi,yj)h2 
ij 


E1  Uu(x0,yj)  -  u(xi,yj)  |  ^  Uu(xn,yj)  -  u(xn-i,yj) 

h  +  2^  1  h 

3= 1  3=1 

,  u(xi,y0)  -  u{xi,yf)  |  u(xi,  yn)  -  u{xu  y 

2^  1  h  +  2^  h 


n—  1 


) 


tan 


du 

dn 


ds  +  O(h). 


On  the  other  hand,  by  integrating  both  sides  of  the  partial  differential  equation  (16)  with  /  =  0  and 
u o  =  0,  we  get 


fl  A  udxdy  =  ff  (  f  v(s)5  (x  —  X(s))  (y  —  Y(s))  dsj  dxdy, 

J  Jn  J  Jn  \J r  J 

or  equivalently, 

[  —ds  —  [  v(s)ds. 

Jan  dn  Jr 

This  completes  the  proof. 

Remark  3.3.  The  double  integral  ffn  Audxdy  can  be  divided  into  three  parts 

/  /  Audxdy  =  /  /  Audxdy  +  Audxdy  +  Audxdy 

J  J Q,  J  J J  J  J Q,e 

du 


du  du  .  „„ 

/  -3T- ds  -  /  T7-ds+  /  TT-ds, 

Ian  on  Jr+  dn  Jr+  dn 
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see  Fig.  1  for  an  ilhisfration,  see  also  [7].  As  e  — >  0,  we  have 

lim  /  /  A udxdy  =  /  v (s)  ds  = 

e^°JJne  Jr  Jr 


du 


ds 


from  the  partial  differential  equation  (16)  with  /  =  0  and  uq  =  0.  Thus  we  get 

du 

Hi,  Jan 


lim  /  /  A  udxdy  =  I  ~—ds. 

.1  In  Ion  071 


Figure  1:  A  diagram  of  the  domain,  interface,  and  integration. 


Lemma  3.4.  There  is  a  solution  to  the  point  source  problem 


AG  (x,xim)  =  <5(a;  -  xi)S{y  -  ym), 


G  (x,xim) 


=  0. 
an 


(30) 


Proof:  We  know  that  iti(x)  =  log  |x— x;m|/(27r)  is  a  solution  to  Aui  (x,  x;m)  =  5{x—xi)8(yy—ym). 
Let  U2  (x,  x;m)  be  the  solution  to  the  following  problem 

Au2  (x,xim)  =  0,  u2  (x,xim)|on  =  Ui  (x,xim)|an  .  (31) 

This  is  a  regular  problem  and  the  solution  u2  is  unique.  Then  G  (x,  x;m)  =  U\  —  u2  is  the  solution 
to  (30).  Note  that  Gh(xij,  x;m)  is  a  discrete  approximation  to  G(x,j  ,  x;m).  The  Green  function 
AG  (x,  x;m)  is  a  little  different  from  the  standard  Green  function  G  (x,j ,  x;m)  for  which  the  boundary 
condition  is  undefined  and  the  solution  can  differ  by  a  constant. 

We  also  define  an  interpolation  function  for  G,l(x,j ,  x;m),  for  example, 

Ghj  (x(s),xim)  =  ^2  h2Gh(xij,xim)Sh  (xi  -  x)Sh  (yj  -  y) .  (32) 

ij 

Note  that  the  coefficients  of  the  interpolation  are  of  0(1).  Other  type  interpolation  schemes  such 
as  bi-linear  interpolation  will  lead  to  the  same  conclusion. 

By  similar  procedure  in  proving  the  discrete  Green’s  formula,  we  can  get  the  second  discrete 
Green’s  formula. 


Lemma  3.5.  Assuming  that  the  distance  between  F  and  90  is  0(1),  that  is,  dist(T,d£l)  ~  0(1), 
then  we  have, 

A hu(xi,yj)Gh(yLij,xirn)h2  =  J  v(s)G(X(s),  x/m)ds  +  0(ft). 

r 


Proof:  Again,  we  show  the  second  discrete  Green’s  formula  by  expanding  the  summation.  After 
cancellation  of  interior  terms,  only  boundary  terms  and  a  source  are  left.  Thus,  we  get 

Y.  &hu(xi,  yj)Gh(xij,xim)h2  =  ^  Gh(xx  j,xim) 

ij  i= 1 

,  ,.«(*«.%•) -“(^n-i,^-)^/,/..  ..  x  ,  V'1  Uu(xi,y0)  -  w(a;*,2/i) ^ 

+  /  v  h  7  ^  v^n_  1  ,j7  5  X/m j  H-  y  ^  tl  ij 


3= 1 
n—  1 


i=l 


,  u{xi,yn)  -  u(xi,yn--i)  h  ^  G/*(x0j,x;m)  -  G?l(xij,xim) 

>  y  h - 7 - G  (Xi,„_l,X;m)  -  2_^  «■ - - - 7 - - - u(xi  ,yj) 


i= 1 


i=i 


\r^  G/t(x„i,xim)  -  G'l(x„_ij,x/m)  (r-4  G,l(xj;o,Xim)  -  G'l(xa,xim) 

-  2^  - - - 1 - - - u(xn-l,yj)  -  2^  - 7 - U[Xi ,  7/i ) 


i=i 


i=l 


Ej  G  Cr  l ,  X/m)  ,  \  /  \ 

h - - - - - u(Si,  yn—i  )+U{Xi,  Urn) 

i= 1 

=  J  [^7^  G(x,  x;m)  -  ^  ds  +  u(xi,ym)  +  0(h). 

This  is  because  \Gh(xitj,xim)  —  G(x,  x;m)|  <  Ch/dist(xij,xim)  for  boundary  points  x.y  when  xim 
is  near  or  on  the  interface  T  and  far  away  from  the  boundary  90,  where  C  is  a  constant,  see  for 
example,  [1,  4,  5].  On  the  other  hand,  by  integrating  both  sides  of  the  partial  differential  equation 
(16)  with  f  =  0  and  uo  =  0,  we  get 


G(x,  xim)Audxdy  =  //  I  /  v(s)6  (x  -  X(s))  (y -Y(s))  ds  )  G(x,xim)dxdy, 
n  J  Jn  \Jr 


or  equivalently, 

J  ^ G(x, X;m)  -  y  \  ds  +  JJ  MAG(x, xlm)dxdy  =  J  v(s)  G(x, xim)ds. 


Note  that 


uAG(x,  xim)dxdy  =  //  u(x,y)5  (x  -  Xi)  5  (y  -  ym)  dxdy  =  u(xi,ym). 
n  J  Jn 


This  completes  the  proof. 
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Lemma  3.6. 


^C^Gh(xy,xim)/i2  =  J  v(s)G(X(s),xim)ds  +  0(/ilogft). 

ij  r 


Proof: 


Nb 

■£c!f  GV„,*ro)i.2  =  ££  VkSh  (xi  -  Xk)6h  (; yj  -  Yk)  AsfeGft(xy,xim)ft2 

ij  ij  k—1 

Nb  __ 

=  ^  % Asfc  ^  (xi  -  a:)  (yj  -  y)  Gh  (x,j,xIm)  ft2 
k=l  ij 


Nb 

=  ^  ^ ^fcAsfcGj  (Xfc,x/m) 

fc= i 


^  ^  ^fcAs^Gj  (Xfc,  x/m)  +  ^  ^  ^fcAs^Gj  (Xfc,  x/m) 

|Xfc-xtm|>VF/i  |Xfc-xIm|<W7i 


[  v{s)G'}  (X(s),xim)ds  +  0(h)+  f  i>(s)G£(X(s),xim)ds  +  0(ft) 

Jr-r'*  Jr'1 


/  u(s)G(X(s),x/m)ds  +  0(ft)  +  /  v(a)G(X(a),x,m)ds  +  0(ft) 
/r— rh  Jrh 


=  /  u(s)G  (X(s),xim)  ds  +  Q(h log  ft). 


In  the  derivation  above,  r£im  is  the  part  of  the  interface  that  intersects  with  the  ball  |X;,.  —  x;m|  < 
IT  ft,  see  the  diagram  Fig.  2.  In  the  last  two  lines  above,  the  second  integral  is  bounded  by  ft,  log  ft 
since  the  integrand  may  have  a  logarithm  singularity  and  the  integral  interval  is  0{h). 


Theorem  3.7.  With  the  assumptions  about  the  discrete  delta  function,  the  IB  is  first  order  conver¬ 
gent  with  a  logarithm  factor  in  the  L°°  norm 


\Eim\  <  Cftlogft,  l,m=  1,2, •••  ,n-  1. 


(33) 


Proof:  Consider  the  error  at  a  grid  point  £)m,  if  x;m  is  close  to  the  interface,  that  is, 
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Figure  2:  A  diagram  used  in  proving  Lemma  3.6.  The  interface  cuts  through  the  ball  centered  at 
xim  with  the  radius  Wh.  Along  the  piece  of  the  interface,  the  singular  integral  needs  to  be  treated 
carefully. 

dist(r,xim )  <  Wh ,  we  have 

Elm  =  {{An)-1  TIB)k[ 

=  ((Ah)~1TIrfg)  kl  +  ((Ah)-1TIi%)  im 
=  0(h2)  +  ((Ah)~1TIi^r)  lm 

=  ^2  fh2Tij(Ah)  1ey—  +0(h 2) 

dis(xy,r)<W/i  '  '  lm 

=  y]  h2  (A  hu(xi,yj)  -  C-f)  Gft(xy-,x/m)  +  0(h2) 

dis(xijtr)<Wh 

=  ^2  h2 Ahu(xi,  yj)Gh(yiij ,  xim)  -  ^  h2C-fGh(xi:j,  x*m)  +  0(/i2) 
ij  ij 

u(s)G(X(s),  x/m)ds  -  ^  ufeG?  (Xfc),  xim)Asfcj  +  0(/i) 

=  0(/ilog/i), 

after  we  apply  Lemma  3.5  and  Lemma  3.6.  Note  that,  in  the  expansion  of  the  summation  from 
dis(xij,T)  <  Wh  to  all  the  interior  grid  points,  we  have  used  the  fact  that  A hu(xi,  yj)  =  0(h 2)  and 
C'2  =  0  when  disfaj,  F)  >  Wh.  If  dist(T,xim)  >  Wh,  the  proof  above  is  still  true  except  that  we 
are  not  going  to  have  a  singular  integration.  Thus,  we  do  not  need  the  log  h  factor. 
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4  Conclusions  and  acknowledgments 


We  give  a  convergence  proof  of  the  immersed  boundary  (IB)  method  in  the  L°°  norm.  We  show 
that  the  IB  method  is  indeed  first  order  accurate  with  a  log  h  factor  if  a  reasonable  discrete  delta 
function  is  used. 
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